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1.  latroductloti 

^This  work  was  undertaken  to  learn  reasons  for  the  behavior  of  the 
^*barrel-stave*^ acoustic  radiator,  operating  in  a condition  devoid  of  all 
pressure  release.  This  source  configuration  is  most  difficult  to  deal 


with  because  it  has  a shape  which  is  mathematically  awkward.  It's 
height  is  small  in  comparison  with  the  diameter,  thus  making  cylindri* 


cal  coordinates  of  little  value.  Admittedly,  the  cylindrical  radiator 


does  not  closely  approximate  a sphere,  but  a circumscribing  sphere 
(radius  the  same  as  the  source  radius)  has  distinct  advantages.  The 
radiation  from  a spherical  source  with  an  arbitrary  surface  velocity 
distribution  is  easily  handled  analytically.  It  is  this  fact  that 
makes  the  present  approach  attractive.  Also,  some  counter  results 
provide  plausible  explanations  for  unexpected  phenomena  that  occur  in 
tests  with  these  radiators. 

In  the  interest  of  completeness,  the  theory  of  the  sj^rlcal  radia- 
tor will  be  presented.  Applications  of  this  theory  to  the  special  problem 
at  hand  indicate  what  can  be  expected  from  the  method.  Consideration  of 
different  frequencies  and  different  velocity  distributions  provide  con- 
clusions that  are  very  useful  in  the  basic  design  process  and  in  the 
process  of  modifying  the  performance  of  existing  designs,  rv 


W 
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11,  The  Theory  of  Radiation  From  a General  Spherical  Source 

The  general  solution  of  a spherical  radiator  can  be  found  in  several 
references » (1),  (2),  (3),  For  the  application  at  hand,  it  will  be  derived 
in  more  detail  to  unify  the  analytic  procedure  and  symbols. 

Let  a spherical  radiator  with  radius,  a,  be  innersed  in  an  infinite. 
Isotropic  and  homogeneous  fluid  of  density  p and  sound  velocity  c.  The 
acoustic  field  outside  the  radiator  can  be  expressed  by  a scalar  wave 
equation  in  spherical  coordinates,  where  t is  the  velocity  potential,  and 
k - u)/c. 


a . ,2  sfi  , _1_  8 Sij 

2 ^^2  a 2 sin  0 W-  ® be-*  , 2.  ^2  * ’ 

c dt  dr  sin  0 ^ 


(1) 


If  there  is  no  velocity  potential  variation  with  the  angle  0 (which 
corresponds  to  the  breathing  radial  mode  of  the  cylindrical  radiator) , 
then  eq.  (1)  can  be  simplified  to  the  two  dimensional  case. 


“2  H " ^ ^ 7ihs  ® 


(2) 


The  solution  of  this  equation  is  well  known  and  can  be  easily  found  to  be: 

-iu)t 


♦ - ) A P (cos  0)C  J (kr)  + in  (kr)]  e" 
Ci  m m m m 

oK) 


(3) 


are  constants  to  be  determined  by  the  boundary  conditions.  The 


quantities  I’lgCcos  0),  Jg,(hr) , and  n.^(hr)  are  Legendre  polynomials. 


spherical  Bessel  functions,  and  spherical  Nevnann  functions  respectively. 

Since  the  radial  component  of  velocity  is  the  radial  derivative  of 
the  velocity  potential: 
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Ur(r,«)  - |i  - 1^1^  - J*.  V'”  9)k[Ji(kr)*l»;(kr)].-‘“‘  (4) 

mH) 

A given  velocity  distribution  on  the  surface  of  the  spherical  radiator 
such  as  U(a»6)e~^'"^  can  be  expanded  in  terras  of  a series  of  Legendre 
functions  of  the  form: 


U(a,0) 


>iu)t 
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with 


^ U P (cos  0) 
•-.ram 

m"0 


r 

U (a,0)  - (m  + %)  u(a,0) 

m Jq, 


P (cos  0)  sin  0 d 0 
tn 


(5) 


(6) 


Where  ra  is  an  integer  and  o and  0 are  the  values  of  0 between  which  the 
active  surface  exists. 

When  r equals  a,  eq.  (4)  should  meet  the  boundary  condition 
of  eq.  (5).  By  equating  them, 


U 

A (a,0) S (7) 

feCj*(ka)  + in»(ka)] 

Q ID 


From  the  results  of  eqs.  (3),  (7),  the  following  acoustic  field  para- 
meters can  be  obtained. 


The  resulting  outgoing  pressure  wave  will  be: 


P ■ P M ■ -impt  - -iiup  y A P„(cos  0)[j  (kr)  +in  (kr)]e"^‘“’^  (8) 

ot  4-  m m m ra 
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and  the  acoustic  impedance  may  be  expressed  as: 
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wtu)p  y A P (cos  e)[j  (kr)  ♦In  (kr)]2iir*8in  6 dee 
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Yap  (cos  e)  kCj'(kr)  +ln*(kr)] 
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(10) 


The  acoustic  Impedance  presented  to  the  radiator  can  be  obtained 
by  setting  r - a in  equation  (10). 

Based  on  the  equations  (4) , (8) , and  (10) , other  acoustic  field 
quantities  such  as  intensity  and  radiated  power  can  be  found.  The 
Intensity  is  given  by  the  real  part  of  the  product  % P*  (where  p^ls 
the  complex  conjugate  of  p).  The  radiated  power  is  determined  by  inte> 
grating  the  intensity  over  a closed  stirface.  At  the  present » however, 
interest  will  be  centered  on  the  development  of  a computer  program  which 
will  Itfc  eentered  on  the  development — ef-e  computer  program  which  will  > 
calculate  the  far  field  pressure  pattern  of  the  short  cylinder. 


III  The  Program 

The  last  equatlon^for  pressure  is  the  one  that  has  been  programmed 
for  computer  solution.  An  effort  was  made  to  make  the  computation  time 
small  so  that  excessive  computer  time  would  not  be  necessary.  It  is  not 
the  ultimate  in  compactness  and  spped  as  it  might  have  been  had  it  be 
deribed  for  a solution  to  one  specific  problem.  It  has  been  deliber- 
ately made  as  general  aB  possible  within  the  confines  of  the  problem  that 
was  initially  specified.  It  contains,  for  example,  numerous  sense  swltcdi 
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settings  enabling  the  user  to  check  the  progress  and  accuracy  of  a 
solution  at  many  points  within  the  body  of  the  main  program.  Instruc- 
tions for  additional  computations  can  easily  be  Inserted. 

The  program  has  restrictions,  of  course,  but  they  do  not  in  any 
way  Impede  Its  intended  use.  All  of  the  limits  and  restrictions,  which 
are  sunnarlzed  later  in  this  report,  can  be  easily  Increased  to  permit 
investigation  of  durface  velocity  distributions  of  greater  complexity 
than  those  thus  far  attempted. 


i 


6 


■ BT” 


Analyst*  of  Results 

I.  In  anelyxing  the  conqiuter  produced  pressure  patterns  it  1*  necessary 
to  keep  in  nlnd  that  a spherical  coordinate  system  was  employed  in  the 
development*  Consequently , the  analysis  will  be  in  terms  of  an  equivalent 
circumscribed  spherical  radiator. 

The  actual  cylindrical  radiator  which  provided  the  experimental  data 
was  approximately  four  feet  in  diameter  and  had  an  active  axial  length  of 
approximately  tticnty  inches.  The  circumscribed  sphere  was  arbitrarily 
established  as  having  a two  foot  radius. 

The  patterns  produced  many  interesting  and  not  altogether  unexpected 
results  from  which  certain  interpretations  and  conclusions  will  be  drawn. 

It  is  to  be  noted  that  in  describing  any  surface  velocity  distribu- 
tion in  terms  of  Legendre  polynomials,  the  closeness  of  "fit"  depends 
largely  upon  the  physical  size  of  the  sphere  relative  to  the  operating 
wavelength.  For  this  purpose  the  parameter  (ka)  is  used  in  which  (a)  is 
the  radius  of  the  sphere. 


The  sphere  is  a geometric  shape  with  many  possible  modes  of  vibration. 
Theae  aiodes,  or  resonances,  determine  to  a large  extent  the  magnitudes 
of  the  various  Legendre  terms.  In  other  words,  a Legendre  term  containing 
a frequency  term  idiose  wavelength  fits  the  critical  dimension  of  the  sphere 
will  be  enhanced  whereas  a term  which  does  not  fit  as  well  will  be  dimin- 
ished. Certain  tersm  may  be  rejected  altogether  if  the  sphere  presents 
to  them  an  anti-resonance.  In  most  cases  this  results  from  the  nature  of 
the  driving  disturbance,  a completely  even  synnstry  in  the  velocity  distri- 
bution giving  rise  to  only  even  order  Legendre  tersm  and  a coaqtletely  odd 
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source  distribution  giving  rise  to  only  odd  terns.  Symnetry  in  this 
sense  refers  to  syooletry  about  the  equator  or  about  the  angular  position 


In  the  following  discussion  a %X  resonance  refers  to  that  condition  of 
operation  which  the  circumference  of  the  sphere  is  one  wavelength,  or  ka"l. 
Similarly,  a ^ X resonance  implies  that  ka~3.  For  the  values  of  ka  less 
than  1,  the  predominant  Legendre  term  is  P^,  ^ich  is  a constant  indepen- 
dent of  the  polar  angle.  This  term  characterizes  the  behavior  of  the  sphere 
(or  other  shapes)  at  extremely  low  frequencies  for  which  ka  < < 1.  Since 
a given  Legendre  polynomial  contains  components  of  all  lower  order  terms 
possessing  the  same  synnietry  no  pattern  is  the  result  of  a single  term  with 
the  exception  of  the  P^  and  P^^  polynomials  which  contain  only  one  term  each. 
Operation  at  exactly  one  of  the  resonance  modes,  however,  does  produce 
patterns  which  are  controlled  essentially  by  a predominant  single  term. 

In  selecting  velocity  distributions  it  is  to  be  noted  that  segments 
of  constant  amplitude  were  employed.  This  was  done  for  several  reasons. 

The  computer  program  is  siaq>lified  in  that  portion  which  evaluates  the 
coefficients,  A^,  from  the  velocity  distribution.  In  addition,  the 
"filtering  effect"  of  the  sphere  resonances  causes  in  general,  the  lower 
order  terms  of  the  Legendre  polynomials  to  predominate  in  determining  the 
far  field  pressure  effects.  In  other  words,  the  idealized  distribution 
shown  in  fig.  1 can  contribute  primarily  those  components  which  coincide 
with  the  first  few  modes  as  determined  by  ka. 

11.  It  can  be  seen  in  fig.  1 that  with  the  velocity  distribution  given  and 
with  ka"l  that  the  dominant  mode  of  operation  is  predetermined  by  the  even 
symmetry  of  the  velocity  function  to  consist  of  even  order  Legendre  poly- 
nomials only.  The  symmetry,  as  mentioned  previously,  is  referred  to  an 
acoustic  axis  taken  at  polar  angle  9 ■ ^ . 
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This  driving  mode  Is  then  a radial  motion  o£  a belt  around  the  sphere 

IT  2ir 

between  the  polar  angles  — and  ~ radians.  It  is  to  be  noted  that  although  1 

the  sphere  is  in  resonance  In  the  mode,  no  output  is  observed  due  to  this  j 

1 ^ i 

mode  since  the  resonance  must  be  excited  by  on  odd  symmetric  velocity  < 

function.  ' 

' i 

The  nearest  acceptable  modes  are  the  and  the  pattern 

I 

due  to  the  P^  would  be  a circle,  and  with  respect  to  a circle  drawn  in  for 

2 

coag>arl8on,  the  effects  of  a slight  j (X)  (two  halves  lambda)  mode  can 

be  seen.  The  coefficients  A , of  the  Legendre  terms  are  determined  princi- 

m 

pally  by  the  (j*(ka)  + i n*(ka))  in  tlie  denominator  of  the  A expression. 

' •'m  m m 

Sines  the  spherical  Bessel  derivatives  increase  rapidly  with  the  order  m for 

I 

small  argument,  ka,  the  higher  order  terc^  in  disappear  equally  rapidly. 

Therefore,  for  the  pattern  under  discussion  the  only  terms  of  significance 

are  the  P and  the  P>  . In  otlier  words,  the  sphere  tends  to  breathe  as  a I 

0 2’*^  j 

unit  and  the  resulting  far  field  pressure  is  quite  uniform  as  a function  | 

I 

of  the  polar  angle.  As  (ka)  is  increased  the  effect,  of  the  next  higher  j 

order  mode,  P^,  is  observed.  Due  to  the  effect  of  the  C j^(kr)+in^(kr)]  | 

term  the  far  field  approximation  causes  each  succeeding  Legendre  term  to 

be  of  alternate  sign,  i.e.  the  terms  add  as  P^  - P2  P^  - P^  etc.  Since 

this  series  of  terms  is  convergent,  however,  each  succeeding  term  generally 

contributes  less  than  its  predecessors  with  the  important  exception  of 

resonant  operation.  In  the  cases  studied  with  (ka  ■)  up  to  4.5  and  with 

either  pure  even  or  pure  odd  excitations  the  pattern  can  be  interpreted 

reasonably  well  by  observing  notch  locations  and  relative  lobe  strengths 

and  comparing  these  with  hand  sketches  of  various  Legendre  combinations. 

The  first  interesting  departure  from  the  pure  breathing  mode  of 
operation,  which  incidentally  implied  rigid  caps,  allows  the  polar  caps 
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to  move  such  that  a constant  volume  is  maintained.  The  pattern  character- 
istically assumes  the  form  of  fig.  6 and  7 with  a sharply  defined  notch. 

Any  pattern  with  deep  notches  cannot  contain  the  term  in  any  appreciable 
SEoplitude.  Since  the  excitation  is  still  of  even  symmetry,  tlie  first  allow- 
able term  is  but  with  ka**!  for  this  figure  there  can  be  little  output 
so  far  from  the  rescnance  mode  (or  ka°>2). 

The  normal  convergence  of  the  coefficients  at  ka"l  precludes  any 

significant  output  at  any  higher  mode,  hence  the  pattern  is  almost  a pure 

P2  but  of  very  low  amplitude.  Increasing  (ka)  to  2.5  shows  essentially 

the  same  content,  but  at  a higher  output  level.  Actually  at  ka"2.5  the 

2 

resonance  point  for  the  mode  has  been  passed.  This  occurred  at  ka~2 
and  would  have  yielded  a maximum  output  if  operated  at  the  corresponding 
frequency. 

The  next  significant  patterns  occur  whenever  an  odd  mode  of 
excitation  is  employed,  see  fig.  8 and  10.  In  this  mode  there  will  always 
be  a complete  notch  along  the  active  or  acoustic  axis.  With  odd  excitation, 
however,  the  Legendre  terms  will  be  restricted  to  the  odd  orders.  The 
lowest  odd  order  is  the  P^  which  would  be  in  resonance  at  ka**l  correspond- 
ing to  the  ^X  mode.  A resonance  pressure  pattern  at  ka**!  would  be  almost 
a dipole  pattern  with  the  major  radiation  occurring  along  the  polar  axis. 

The  actual  barrel  stave  transducer  to  which  reference  has  been  made  previous- 
ly is  built  in  two  identical  sections  so  that  there  are  in  reality  two 
active  bands  of  barrel  staves  adjacent  to  each  other.  The  results  of  the 
pattern  of  fig.  8 indicate  that  instead  of  driving  the  two  section  in  phase, 
as  has  been  customary,  opposite  phased  excitation  would  give  rise  to  a strong 
dipole  pattern  vertically  oriented.  The  actual  (ka)  used  in  the  patterns 
of  fig.  8 was  2.5,  hence  the  appearance  of  a small 
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dipole  pattern  is  dravm  for  comparison,  and  it  can  be  seen  that  at  ka~2.5 
the  output  is  essentially  "dipole  like". 

The  excitation  function  for  fig.  10  is  for  a double  reversal  of 
velocity  along  the  active  length  and  the  "filtering"  action  of  the  sphere 
maintains  essentially  the  same  Legendre  content  until  kas'S  where  further 
notching  will  appear. 

It  is  postulated  here  that  the  velocity  distribution  of  figures  8,  9, 
and  10  are  possible  if,  at  certain  critical  frequencies,  a flexural  standing 
wave  exists  along  the  length  of  the  stave.  If  so,  this  would  give  rise  to 
the  patterns  of  the  above  mentioned  figures  depending  upon  the  boundary 
conditions.  The  velocity  distribution  of  figure  9 is  that  of  figure  8 phase 
shifted  along  the  length  of  the  stave.  The  excitation  for  the  pattern  of 
figure  9 is  therefore  even  symmetric  . of  the  type  principally,  but  with 
a reduced  output  along  the  acoustic  axis.  This  is  comparable  to  fig.  12 
redrawn  from  test  data.  This  obviously  is  a primary  excitation  with 
some  P^  superimposed. 

Particular  attention  should  be  drawn  to  the  pattern  of  figure  9, 

%(hlch  suggests  that  a flexural  wave  might  be  the  cause  of  sudden  reductions 
in  acoustic  axis  output  at  a particular  frequency.  This  form  of  flexural 
wave  could  exist  if  the  flexural  wave  velocity  were  in  the  neighborhood 
of  2000  ft/sec.  Operation  at  a (ka)  value  of  2 would  seem  to  produce  a 
reasonable  expectation  of  this  mode  of  operation. 

Preliminary  tests  on  a single  stave  indicate  a flexural  wave  velocity 
in  the  right  order  of  magnitude.  At  the  moment,  however,  the  prime  concern 
la  the  interpretation  of  pressure  patterns  on  the  basis  of  assumed  velocity 
distributions  and  not  on  the  detailed  reasoning  for  the  existence  of  such 
distributions. 
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One  ccoibination  o£  breathing  mode  and  a superimposed  flexural  wave 
has  been  included  in  this  report,  see  fig.  11.  It  indicates  a simple 
addition  of  the  individual  patterns.  Other  combinations  immediately  come 
to  mind  and  some  of  these  are  currently  under  study  and  will  be  reported 
at  a later  date. 
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Conclusions 

1.  The  spherical  approach  to  the  cylindrical  radiator  seems  reasonable 
in  so  far  as  it  has  been  applied  to  specific  pressure  patterns  of 
an  axially  short  cylinder. 

2.  To  produce  a breathing  mode  of  any  significant  value  requires  that  the 
real  or  assumed  velocity  distribution  have  an  average  value  over  the 
spherical  surface. 

3.  Spherical  resonances  appear  to  control  the  pattern  shape  which,  for  the 
usual  range  of  ka  values,  are  dictated  by  a relatively  few  low  order 
modes. 

4.  The  existence  of  a flexural  mode  will  explain  certain  pattern  notches 
or  other  sudden  reductions  in  acoustic  output  along  a particular  axis. 

5.  Multi>section  transducers  could  be  phased  to  produce  changes  in  direcv 
tivity  such  as  the  dipole  effect  discussed  previously. 

6.  Probable  surface  velocities  can  be  ascertained  by  observation  of  the 
far  field  pressure  pattern. 

7.  The  influence  of  volume  effects  on  the  pressure  pattern  can  be 
observed. 

Many  workers  in  the  field  of  design  of  large  radiators  have  used 
similar  mathematical  approaches  to  a model  of  a short  pressure 
equalized  cylinder.  Some  of  these  seem  much  more  appropriate  than 
does  a circumscribed  sphere.  One  such  is  the  prolate  spheroidal 
coordinate  system.  The  real  Justification  for  the  spherical  model  is 
two-fold,  however.  The  associated  mathematics  is  vastly  simpler  on 
the  one  hand,  and  there  are  only  very  minor  differences  between  the 
corresponding  source  distributions.  That  is,  the  sphere  and  the 
prolate  spheroid  are  so  similar  and  so  close  to  each  other,  that  only 
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minor  differences  in  particle  velocity  are  possible. 

Clearly,  these  advanteges  of  the  sphere  may  be  lost  when  the  height 
to  diameter  ratio  excedes  unity.  But  when  it  is  appropriate,  the 
spherical  model  makes  possible  an  appreciation  of  the  individual  contri- 
butions of  various  structural  modes  to  the  net  radiation  and  radiation 
impedance  effects.  Knowing  ehat  the  individual  modal  contributions  are 
makes  it  possible  to  enhance  or  suppress  certain  ones  so  as  to  modify 
the  acoustic  properties  of  a specific  design. 

Finally,  it  may  be  said  that  the  present  design  approach  is  an 
over-simplification  which  is  warranted  because  of  the  Insight  it 
provides,  even  though  the  quantitative  results  are  inherently  in- 
accurate. The  cooqputer  results  so  far  have  provided  logical  reasons 
for  the  appearance  of  notches  in  the  frequency  response  curves  and 
the  effect  of  rigid  sections  in  the  polar  cap  regions  (top  and  bottom) 
of  the  cylindrical  transducer. 
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